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Abstract 

We examine the two elliptic systems given by 

{G) Xn ~Au = Xf'(u)g(v), -Av = jf(u)g'(v) in O, 

and 

(H) Xn -Au = \f(u)g'(v), -Av = r f(u)g(v) in 

with zero Dirichlet boundary conditions and where A, 7 are positive parameters. We show that for 
arbitrary nonlinearities / and g that the extremal solutions associated with (G)a, 7 are bounded provided 
Q is a convex domain in where N < 3. In the case of a radial domain we show the extremal solutions 
are bounded provided N < 10. 

The extremal solutions associated with (H)\ tl are bounded in the case where / is arbitrary, g(v) = 
(v + l) 9 where 1 < q < 00 and where SI is a bounded convex domain in R N , N < 3. 

Results are also obtained in higher dimensions for (G)a, 7 and (H)\ tJ for the case of explicit nonlin- 
earities of the form f(u) = (it + l) p and g(v) = (v + l) q . 

2010 Mathematics Subject Classification. 
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1 Introduction 

In this work we examine the following systems: 



C -Au. = \f{u)g{v) il 
(G)a, 7 <^ -Av = jf(u)g'(v) 0, 

[ u = v = on 

and 

*This work is supported by a University Graduate Fellowship and is part of the second author's Ph.D. dissertation in 
preparation under the supervision of Professor N. Ghoussoub. 
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-Ait = Xf(u)g'(v) Q 
(H) x „ { -Av = if{u)g{v) n, 
u = v = <9f2 

where f2 is a bounded domain in M. N and A, 7 > are positive parameters. The nonlinearities / and g will 
satisfy various properties but will always at least satisfy 

(i?) / is smooth, increasing and convex with f(0) — 1 and / superlinear at 00. 

We begin by recalling the scalar analog of the above systems. 
Given a nonlinearity / which satisfies (R), the following equation 



(Q)> 



-An = A/O) n 
u = an 



is now quite well understood whenever f2 is a bounded smooth domain in Mr. See, for instance, [2j|3l|4l[9j 
HU [121 [1] . We now list the properties one comes to expect when studying (Q)\. It is well known that 
there exists a critical parameter A* € (0, 00), called the extremal parameter, such that for all < A < A* 
there exists a smooth, minimal solution u\ of (Q)\. Here minimal solution means in the pointwise sense. In 
addition for each x G the map A t— > u\(x) is increasing in (0, A*). This allows one to define the pointwise 
limit u*(x) := lim,\/<A* u\(x) which can be shown to be a weak solution, in a suitably defined sense, of 
(Q)x' ■ F° r this reason u* is called the extremal solution. It is also known that for A > A* there are no weak 
solutions of (Q)\. Also one can show the minimal solution ux is a semi-stable solution of (Q)\ in the sense 
that 

/ a/'M</> 2 < f iwi 2 , y^eH^n). 

Jn Jn 

A question that has attracted a lot of attention is the regularity of the extremal solution. It is known that 
the extremal solution can be a classical solution or it can be a singular weak solution. We now list some 
results in this direction: 

• u* is bounded if f(u) = e u and N < 9. 

• (El) u* is bounded if / satisfies (R) and N < 3. 

• ([3]) u* is bounded if / satisfies (R) (can drop the convexity assumption) and a convex domain in 
R 4 . 

• ([!]) u* is bounded if fl is a radial domain in M. with N < 10 and / satisfies (R) (can drop the 
convexity assumption) . 

It is precisely these type of results which we are interested in extending to the case of systems. Before we 
can discuss the regularity of the extremal solutions associated with (G)a, 7 and (H)\ n we need to introduce 
some notation. 

Under various conditions on / and g the above systems fit into the general framework of developed in 
[13) , who examined a generalization of 

( -Ait = XF(u,v) 

[ u = v = on. 

The following results are all taken from [T3]. Let Q = {(A, 7) : A, 7 > 0} and we define 
U := {(A, 7) G Q : there exists a smooth solution (u,v) of (-P)a,7} • 
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Firstly we assume that F(0, 0), G(0, 0) > 0. A simple argument shows that if F is superlincar at u = oo, 
uniformly in v, then the set of A in hi is bounded. Similarly we assume that G is superlincar at v — oo, 
uniformly in u and hence we get hi is bounded. We also assume that F, G are increasing in each variable. This 
allows the use of a sub/supersolution approach and one easily sees that if (A, 7) € hi then so is (0, A] x (0, 7]. 
One also sees that hi is nonempty. 

We now define T := dhi D Q, which plays the role of the extremal parameter A* . Various properties of T 
are known, see [13]. Given (A*, 7*) £ T set a := jj G (0,oo) and define 

T CT :={(A,A(t) : y < A< A*}. 

We let (u\,v\) denote the minimal solution (P)x,a\ for ^- < A < A*. One easily sees that for each x e ft 
that u\(x), v\{x) are increasing in A and hence we define 

u*{x) := lim u\(x), v*(x) :— lim v\(x), 
xy\' a y\* 

and we call (u*,v*) the extremal solution associated with (A*, 7*) G T. Under some very minor growth 
assumptions on F and G one can show that (u*, v*) is a weak solution of (P)>* ~*. 
We now come to the issue of stability. 

Theorem A. Let (A, 7) € hi and let (u,v) denote the minimal solution of (P)a i7 . Then (u,v) is 
semi-stable in the sense that there is some smooth < (,x 6 -ffg(Q) and < r\ such that 

- AC = \F u (u,v)C + \F v (u,v) X + ri(, -A x = 7G u (u,v)C + jG v (u,v) X + VX, (1) 

In this paper we prove that the extremal solution of (G)a*, 7 * with general nonlinearities, either on a 
general domain and lower dimensions or on a radial domain and higher dimensions are regular. Moreover, 
for explicit nonlinearites we prove regularity on a general domain in higher dimensions. 

The following stability inequalities play a key role in this paper and we shall refer to them many times 
through proofs. We mention that in |10j the De Giorgi type results and Liouville theorems have been proved 
for a much more general gradient system and they obtained a stability inequality which reduces to ([2]) in 
the particular case we are examining. Note that some of our results will hold for a general gradient system 
that is examined in |10) . 

Lemma 1. 1. Let (u,v) denote a semi-stable solution of (G)a, 7 in the sense of (QJ). Then 

f"{u)g{v)tf + J f(u)g"{v)^ 2 + 2 J f\u)g'(v)H < ± j |V0| 2 + - j (2) 

2. Let (u,v) denote a semi-stable solution of (H)\^ in the sense of (QJ). Then 

1 /■,_,„. 1 



f'(u)g'{v){4> 2 + V 2 ) + 2 j Vff Jr gg 1J H' <jj |V<^| 2 + -J |V^| 2 . (3) 

Proof. (1) By Theorem A there is some < C, X such that 

-AC > Xf"(u)g(v)( + \f'(u)g'(v)x and -A X > 7/ W(»)C + 7/(X)<?»X in fi, 

Consider test functions 4>,ip E Hq(£1) and multiply both sides of the above inequalities with and — to 

1 s X 

obtain 

|V0| 2 ^ + 2 |V0.VC^>| Xf(u)g'(v)^+ J Xf(u)g"(v)ct> 2 , 

|W| 2 ^ + 2 / • V X 4 > / if(u)g"(vH> 2 + / 7/WW 2 -. 



ip J V 1 
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Apply Young's inequality for the left hand side of each inequality and add them to get 

A J f"{u)g{v)<p 2 + 1 J f(u)g"(v)iP 2 + J f(u)g'(v) (V| +7 ^) < J |V^| 2 + J |Wf, 
Simple calculations show that the third term is an upper bound for 

2^ / f'{u)g'{v)U- 



Then, replacing (f> with and ip with gives the desired result. 

(2) Proof is quiet similar to (1). By Theorem A there is some < £, \ such that 

> \f(u)g'(v) + \f(u)g"(v)~ and - — > jf"(u)g(v)~ + jf'(u)g'(v) in £1, 

C C x x 

and we now multiply the first equation by 4> 2 and the second by ^ 2 and add the equations and integrate 
over fi. In addition we use the fact that 

AE - 2 < w, 



for any E > and <fi £ Doing this one obtains 



f(u)g'(v)(Xcb 2 +^ 2 ) + I Xf{u)g"{v)^ +1 f"{u)g{v)^< f |V0| 2 + |V^| ; 

C X Jn 



in Jn 
Again some simple algebra shows that 



2^ / y/f(u)f»(u)g{v)g»(v)W, 



is a lower bound for the second integral. Using this lower bound and replacing <\> with and ip with 
finishes the proof. 

□ 

In Section 2, we explore the regularity of extremal solutions for systems (G)a, 7 and (H)\.y with arbitrary 
nonlinearities and, in then Section 3 we consider explicit nonlinearities. We finish the current section by this 
point that in [5] the system 

-Ait = Xe" fl 
(E) Xn { -Av = je u it, 

u = v = <9f2, 

was examined. It was shown that if f2 is a bounded domain in M. where N < 9, then the extremal solution 
(u*,v*) associated with (A*, 7*) € T is bounded if 

N-2 7* 8 

^ < T- < 



A* N-2 

Note that as one gets closer to the diagonal parameter range 7 = A that better regularity results are obtained. 
At the diagonal the system can be shown to reduce to the scalar equation — Am = \e u . This phenomena will 
also be present in Section 3 where we consider explicit nonlinearities. 
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2 Arbitrary nonlinear it ies 

We begin by examining (G)a, 7 in the case of arbitrary nonlinearities and we show the extremal solutions are 
bounded in low dimensions and our methods of proof are close to |14) . 

Theorem 1. Suppose that Q, is a bounded smooth convex domain in M. N where N < 3 and suppose f and g 
both satisfy condition (R). We also assume that a :— /'(0) > and b :— <?'(0) > 0. In addition we assume 
that f',g' are convex and there is some £ > such that 



liminf > 0, liminf > 0. (4) 



Let (A*, 7*) G T. Then the associated extremal solution o/(G)a*, 7 * is bounded. 

For radial domains we obtain similar results but in higher dimensions and our methods of proof follow 
very closely to [I] and |15j . 

Theorem 2. Let fi = B\, N > 3, and f and g both satisfy condition (R) and in addition we assume that 
there is some £ > such that 

lim inf > 0, lim inf > 0. 

Let (A*, 7*) G T and let (u*,v*) denote the extremal solution associated with (G),\*, 7 *. Then 

1. if N < 10, then u* ,v* G L°°{B 1 ), 

2. ifN= 10, then u*(r),v*(r) < C A ., 7 »(1 + |logr|) for r G (0,1], 

3. ifN> 10, then u*(r),v*(r) < C x , , 7 » ^r'f +^ 77 ^ T + 2 for r G (0,1]. 

We are unable to prove the analogous version for the system (H)\ t -y and hence we restrict our attention 
to the special case. 

Theorem 3. Suppose il a bounded smooth convex domain in R 3 and 1 < q < oo. Assume f satisfies (R) 
and we also assume that f" > C > 0. Let (A*, 7*) G T. Then the associated extremal solution of (H)\*~ ( * 
for g(v) — (1 + v) q is bounded. 

The following lemma is used to prove Theorem [1] where a convex domain is assumed but we prove the 
lemma for general domains. 

Lemma 2. Suppose ft is a bounded domain in M. N and f and g satisfy the conditions from Theorem\l\ but 
one can weaken condition to f",g" — > 00 at 00. Let (A*, 7*) G T and let (u*,v*) denote the extremal 
solution associated with (G)a*. 7 *. Then there is some C < 00 such that 

(0 J ' f'(u*)9'(v*)(f(u*)-a)(g'(v*)-b)<C, (ii) J (f'(u*) - a)f"(u*)g(v*) < C, 

(Hi) [ (g'(v*) - b)g"(v*)f(u*) < C. 



lim inf ^ > 0) lim inf ^ > Q (5) 



Remark 1. Let fi,gi denote smooth increasing nonlinearities with /i(0), gi(0) > and we also assume there 
is some £ > such 

A^>0, liminf^ 

Let (u m ,v m ) denote a sequence of smooth solutions of (P)\ m , a \ m , where < a < 00 is fixed and X m is 
restricted to a compact subset of (0, 00) and F(u,v) = fi(u)g\(u), G(u,v) — f2{u)g2(v). Then we have the 
estimate 

fi(u m )gi(v m )S + / f2{u m )g 2 {v m )o~ < G, 
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where 8{x) :— dist(x,dQ). Applying regularity theory shows that u m ,v m are bounded in L 1 (f2). 

On occasion we will restrict our attention to smooth convex domains where many of the proofs are much 
more compact. For this we will use a result proven using the Moving Plane Method. So we assume that 
is a smooth convex domain. Then there is some Sq > small such that any smooth solution (u, v) of 
(P)a,7 satisfies some monotonicity properties in in the small strip {x £ : S(x) < Eq} near the boundary. 
(Essentially one can assume u,v are strictly decreasing as one approaches the boundary). This coupled with 
the L 1 bounds on u m and v m shows that there is some constant Co such that u m , v m < Cq on some small 
strip, say {x € fi : 5(x) < £1} near the boundary. Assuming the same restrictions on A m ,7 m one can use 
the maximum principle to show that there is some C\ > such that u m ,v m > C\ in the compliment of this 
small strip. 

Proof. All integrals are over unless otherwise stated. Our approach will be to obtain uniform estimates for 
any minimal solution (u, v) of (G)\ n on the ray r CT and then one sends A /~ A* to obtain the same estimate 
for (u* , v*). Let (it, v) denote a smooth minimal solution of (G)a, 7 on the ray T a and put 4> :— f'(u) — a and 
■ijj := g'(v) — b into ([2]) to obtain 

f"(u)g(v)(f(u) a) 2 + J f(u)g"(v)(g'(v) - bf + 2 J f'(u)g'(v)(f'(u) - a)(g'(v) b), 
is less than or equal 

1 r V(f(u)-a)f"(u)-Vu+- [v(g'(v)-b)g"(v)-Vv. (6) 



Integrating ([5]) by parts shows that this is equal to 

^ J (f(u) a)r(u)\\7u\ 2 + l -J (f'(u) - o)/»(-A«), 

plus a similar term in involving v. We use the equation (G)a, 7 to replace —Au and —Aw in the last line and 
simplify to arrive at 

i J (f(u) - a)/"'( M )|V W | 2 + -J (g'(v) b)g'"{v)\Vv\ 2 + 2 J f(u)g'(v)(f(u) ~ a)(g'(v) - b), 
is less than or equal to 

a { f(u)(f(u)-a)g(v) + b [ g"(v)(g'(v) - b)f(u). 



We now define h x {u) := J^(f'(t) ~ a)f"(t)dt and h 2 (v) := ^{g'{t) - b)g"'{t)dt. Subbing this into the 
previous inequality and integrating by parts and using (G)\ n again we arrive at 



hi(u)f'(u)g(v) + J h 2 {v)f{u)g'{v) 
+2 J f(u)g'(v)(f(u) - a){g'{v) - b) < a J f'(u)(f(u) - a)g(v) 

+b j g"{v){g'(v)-b)f{u) (7) 

Now suppose u > a > 0. Then we have 

hi(u)> / (/' (t) - a)f"{t)dt > (f'(a) - a){f" (u) - f'(a)), 

J a 

and so using the condition on f"(u) we see that 

hminf — — - > / (a) - a, 

«->oo J [Uj 



G 



for any a > 0. But since / is convex and superlinear at infinity we see that limu-^oo jrrj~j = oo. Similarly 
lim^^oo -ttV4 = oo. 



9"(v) 

We now estimate the integral J f"(u)g(v)(f'(u) — a). There is some T > 1 large such that for all it > T 
we have fti(u) > 100(a + l)f"(u) for all u>T. Then we have 



A>t y«<T 100(a+l)7 A<t7 
We now estimate this last integral. Let T be as above and fixed and we let k > 1 denote a natural number. 



f"(u)g(v)(f'(u) -a)= + = C(k, T) + / f"(u)g(v)(f'(u) - a) 

u<T Ju<T,v<kT Ju<T,v>kT Ju<T,v>kT 

and we now estimate this last integral. One easily sees that this last integral is bounded above by 



u<T 



/'(«) :>Tt - %'(«) 



Combining this all together we see that for all sufficiently large T and all 1 < k there is some constant 
C(k,T) such that 

f{u)g{v){f{u)-a) < m{ l +1) J hMf'MgW + CfaT) 

+ sup^M sup [(f'(u)-a)(g'(v)-b)f( U )g'(v). 
u<T f'(u) v> kT (g» - O)ff'(u) J 

Using the same argument one can show for all sufficiently large T and for all 1 < k there is some C(k,T) 
such that 

9"(v)f(u)(g'(v)-b) < —J—j h2 {v)g'{v)f{u) + C{k,T) 

sup ^ SU P 7777 / {U \ flf ^ !{f{u)-a){g'{v)~b)f'{u)g'{v). 



v<T g'{v) u >kT (/'(«) - a)f'(u) 
Since f",g" — > oo we see that 

lim sup M = 0, 

k^oo u>kT {f'{u) - a)f'(u) 

and similarly for the other term. Hence by taking k sufficiently large we can substitute everything back into 
Q and see that all the integrals in are bounded independant of A. 

□ 

Proof of Theorem [TJ 

We assume that N = 3 and f2 is convex domain in R 3 . The case of N = 1 , 2 is easier and we omit 
their proofs. We suppose that (A*, 7*) £ T and (u*,v*) is the associated extremal solution of (G)a* j7 *. Set 
a = Using Remark [T] along with Lemma [5] we see that f'(u*)g'(v*) £ L 2 (f2). Note that this and the 
convexity of g show that 

/Wf < r 

n («* + l) 2 " ' 

From Lemma [5]and Remark [T] we have — Am*,— Av* £ L 1 and hence we have u*,v* £ L 3 - , ie. L p for any 
p < 3. 
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We now use the domain decomposition method as in [TJ]. Set 

O, := fx : > f'(u'f-" s (v'f-* 



y (W*+1) 2 

fi 2 := n\ni = {x : f'(u*)g(v*) < (v* + , 



{f{u*)g{v*)f- a < I f{ f} 29[ X/ < C 



where < a is to be picked later. First note that 

: ft 

(v* + 1) 
Similarly we have 

r (f(u*)g(v*)7< f («* + !)- 



Taking a = | and using the bound on v shows that /' (u*)g(v*) G L 5 -(il). By a symmetry argument 

6 

we also have f(u*)g'(v*) G L 5 -(f2). 

By elliptic regularity we have u*,v* G W ' 5 - and this is contained in L°- (0) after considering the 
Sobolev imbedding theorem. Using these estimates and again using the domain decomposition fii and 2 

but taking a = | gives that f (u*)g(v*) G £ 5 -(f2) and by symmetry we have the same for f(u*)g'(v*). 

2 - 

Elliptic regularity now shows that G TU ' 2 - and this is contained in LP for any p < 00. One last 

iteration with a = k shows that /' (u*)g(v*) G Li(fl) and after considering elliptic regularity and the 
Sobolev imbedding we have u* is bounded. By symmetry we see v* is also bounded. 

□ 

Proof of Theorem [3j Let (A*, 7*) G T and let (u*,v*) denote the extremal solution associated with 
(iJ)A*, 7 *. We let (u,v) denote a minimal solution on the ray F CT where a — Without loss of generality 
we can assume A = a = 1 to simplify the calculations. Note the assumption on /" shows there is some 
£ > such that ^ttt -> 00 as it —> 00 . Using this along with the fact that q > 1 and Remark [T] shows that 
Au*,Av* G L x (d). 

Set a :— f(u) — 1 and (3 := (v + 1)* — 1 where 1 < t < t + (q) := q + J q(q — 1) into the stability inequality 
to arrive at 

te-^rr) / f(u)(v+i) 2t+ o- 1 + j \f(uyi)r(u)\vu\ 2 +2^qJq—T) J v/M/"H(«+i) 9 - 1 (/(«)-i)((^+i) t -i) 

is less than or equal to 

2q J f(u)f(u)(v + iy- 1 +2q J f{u){v + iy +t ~\ 

We label these integrals as Ii for 1 < i < 5 from left to right. The condition on t ensures the coefficient in 
front I\ is positive. We can rewrite 

qhx{u)f(u)(v+iy-\ where hx{u) = /"(/(r) - l)/"(r)dr. 



One easily sees that ^tt~j — ► 00 as u — > 00. Let T be sufficiently large such that hi(u) > 10 f'(u) for all 
u >T. Then one easily sees that 

f(u)f'(u)(v + l)"- 1 < 1 J + l)"- 1 + /'(T) | /(«)(« + I)'" 1 . 

We also have 

/'(«)(« + l)^*- 1 < T^- 1 f /'(«) + -L /" /'(„)(« + 
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and so after combining the estimates we have 

\ ~ - J f(u)(v + l)*-*-x + /*(«)/(«)(« + + 



+2^~ T) y + - 1)((« + 1)* - 1), 

is less than or equal 

2qf(T) [ f(u)(v + iy- 1 +2qT<i+ t - 1 [ /'(«). 



Passing to the limit shows this inequality holds with (it*, v*) in place of (u, v). 

But these last two integrals are finite and hence we have an estimate provided the first coefficient is 
positive, which is indeed the case provided we take T bigger if necessary. Hence each of the following 
integrals is finite 

(i) [ f'(u*)(v* + l) q+2t -\ (ii) f f(u*)f(u*){v* +iy-\ (in) [ f(u*)Hv* + iy +t -\ 



Then we have 



for all 1 < t < t + (q). 

Now note that —A(f(u)) = —f"(u)\\7u\ 2 +f'(u)f(u)(v+l) q ^ 1 and since / is convex and since f'(u)f(u)(v+ 

N 

is uniformly bounded in L 1 (fi) along the ray shows that f(u) is uniformly bounded in L N ~ 2 - = L 3 ~ ; 
and hence f(u*) £ L 3 - . 

We now use (i) and (iii) to show that u* is bounded. Pick p > | but close and a > and 1 < r < oo 
such that 

3 

(p-a)r<3, T 'a = -, (q - l)pr' < q - 1 + t + (q). 

(f{u*)(v* +i) (q - i) ) p < (j j{u*)^- a ^y f{u*)i{v+iy , ^ q - i ^y T > 

but the right hand side is finite and hence by elliptic regularity we have u* is bounded. We now show that 
v* is bounded. First note that we have J(v* + i^t+Q- 1 < oo for any 1 < t < t+(q) and hence one has 
J(v* + l) q+1 < oo. Since u* is bounded this shows that v* £ Hq. Now to complete the proof of v* being 
bounded it is sufficient to show (since u* is bounded) that (v* + 1) £ L^ _1 ^ p for some p > | but this easily 
follows after considering the above estimate. 

□ 

Proof of Theorem [2l Step 1. Let (u,v) denote a smooth minimal solution of (G)\ n on the ray T a 
where a := -pr. Then taking a derivative of (G)a, 7 with respect to r gives 

-Au r + ^u r = \f"(u)g{v)u r + \f(u)g'{v)v r for < r < 1, 
-Av r + ^v r = 1 f l {u)g'{v)u r + 1 f{u)g"{v)v r for < r < 1. W 

Multiply the first and the second equations of (JSJ) with u r (j) 2 and v r <p 2 where <j> £ C 0,1 (i?i) n Hq(B\) gives 

| VUr | 2 2 + U/u 2 r ■ 4? + ^^-ul4> 2 = j \f" (u)g{v)u 2 4 2 + \f'(u)g'{v)v r u r 4> 2 
\Vv r \ 2 4> 2 + \vv 2 r ■ + ^-v 2 r 4> 2 = f 1 f(u)g'{v)u r v r 4> 2 + 1 f(u)g"(v)v 2 ^ 2 (9) 



2 

On the other hand, testing @ on <fi — > u r (f> and ip — > v r <p where <fi is as above, we get 

J f"{u)g{v)u 2 r 4, 2 + J f(u)g"(v)v 2 r <t> 2 + 2 J f(u)g'(v)u r v r ^ < I J |VK0)| 2 + - / V(r,,>,|- 
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Expanding the right hand side we get 

\J (\Vu r \ 2 <f> 2 +ul\V<i>\ 2 + ±V<P 2 -Vui)+l-J (\Vv r \ 2 <j> 2 +v 2 r \Vcf>\ 2 + ±V<t> 2 -Vv 2 r 



A 

Applying © the above will be 



n-i r 2 4> 2 N-i r 2 4> 2 



+2 [ f{u)g'{v)u r v r <p 2 + [ f{u)g"(v)v 2 r tf+ I f{u)g(v)u 2 r tf 



Therefore one obtains, after substituting rip for 



{N ~ 1] I + f ) ^ " / (f + f) |VM|2 ' (10) 
for all G C 0,1 (i?i) n Hq(Bi). Note that there is no / and g in this estimate. 

Step 2. We show that (TlT)|) implies that u*,v* £ Hq(B{). Firstly we argue there is some Cr > such that 
for any < R < 1 we have sup r>R (u(r) + v(r)) < Cr and Cr is independent of A (and hence the estimate 
also holds for u*,v*). To see this we first note that by Remark[T]we have \\v\\li < C (uniformly in A) 

and since u, v are radially decreasing we have the desired result otherwise we couldn't have the L 1 bound. 
We now let < cj) < 1 be a smooth function supported in B\ with <f> = 1 on Bi . Putting this into (flQ|) and 
rearranging gives 



(N-2) ^ + ^< C + (11) 

Jb 1 A 7 Jbi\Bi x 7 

2 2 

Now let < ip < 1 denote smooth function with if) = in Bi and ^ = 1 for |sg| > |. Multiply —Au = 
Xf'(u)g(v) by uif) 2 and integrate by parts and use Young's inequality to arrive at 

|Vw|V < 2A / f'(u)g(v)uif) 2 +4 / w 2 |V^| 2 , 



and hence we have 

/ \Vu\ 2 <2\[ f'(u)g(v)u + C [ u 2 (12) 

JBx\Bi JBi\Bi 

and we now use the pointwise bound to see that 

|Vw| 2 < C, 

Bx\Bi 

2 

where C is independant of A. Similarly we obtain the same estimate of v and combining this with (|llj) we 
see that u, v are bounded in Hq(Bi) indapendant of A and hence u*,v* £ Hq(B±). 

Step 3. Let (u, v) denote a minimal solution of (G)a, 7 on the T a where a :— 
For < r < i define (f> to be the following test function 

{ r -VW=T-i if o < t < r, 

t -VTT=T-i ifr<t<l/2, 

2^=T+ 2 (1 -t) if 1/2 < t < 1. 

Putting cf) into (fTU)) gives 
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for all < r < \ and one easily extends this to all < r < 1 by taking Cn bigger if necesary. From this, by 
simple calculations we get 

Kr)-w(-)| + — |«(r)-«(-)| < / ( — K(t)| + — |« r (t)| ) dt 

1/2 / \ 1/2 



< c^ 2 -f (-L|| Vu || i2(BiTO + -L||v«|| i2(BiTO ) . 

Let < r < 1. Then, there exist m £ N and 1/2 < n < 1 such that r = 2T ^Li . Since u,w are radial, we 
have u(n) < ||«Hloo (ba b^ ^ c, Jv||u|| ff i (BAB ^ and w(n) < IMIz,oo (BABv ^ ^ ^ | M I B i (B ABv 7) ■ 

^N r )l + ^=K r )l ^ -^Hn) - u(r)\ + ^=\v(n) - v(r)\ + ^=\u( n )\ + ^=\v( ri )\ 

^ m— 1 _^ m— 1 

^ 7x§K2&)- w S)l + ^gh(2&)-HS) 

Cat I, || Cat , , 

VA (SASl/2) Vf I|U|Ih1 (Bi\Si /2 ) 

^ / n \ -AT/2+v / iV=T+2 / i i \ 

Cat Cat 

VA (SASl/2) ^7 I|u|Ih1 (Bi\Si/ 2 ) 



/" 1 / n \ -A r /2+ V 0V=T+2 \ / 1 1 



i=i 



Note that the sign of — 1 + 2 — -j is crucial in getting estimates. Since -y/AT — 1 + 2 — y = if 
and only if AT = 10, this dimension is the critical dimension. From the above, for any < r < 1 we get if 
2 < N < 10, 

_L Kr) | + _L Kr) | < Cjv ^||«|| ffl(BA5I ^ + ^IMhhb^—^) > 

if AT = 10, 

+ ^=Kr)| < Cjv(1 + |Iogr|) (^||«|| ff i (BlX B^ + ^||w|Ihi( Bi \b^) , 

if AT > 10, 

Passing to limits we obtain the desired estimates for u*,v* . 

□ 
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3 Explicit nonlinearities 

We now examine the case of polynomial nonlinearities and for this we recall the definition t + (p) = p + 
\/ p(p — 1) and note that t+ is increasing on [1, oo). We begin with the gradient system. 

Theorem 4. Let f(u) = (u + l) p , g(v) = (v + l) q and suppose p, q > 2. 

1. Let (A*, 7*) G T. Then the associated extremal solution o/(G)a*, 7 * is bounded provided 

N 2 

— <1 + — -nuuc{t + (p-l),t + (g-l)}. (13) 

2 p + q - 2 

2. Let{X*,j*) G T and define / p ,,,a, 7 (*) := P + <?~ 1- 5^1 - x ) ■ ie ^o :=max{i+(p- 
l),t + (g — 1)} and suppose that Ip^x* {to) > and L qprl *,\* (to) > 0. 

The map t i-> min{Ip )g) A*,-y*(i), J g ,p, 7 *,A*(£)} *s decreasing and has a root in (to, 00), which we denote 
by T. Suppose 

N 2 

— < 1 + T. (14) 

2 P+q-2 V ; 

T7ien the associated extremal solution o/(G)a*, 7 * is bounded. 

Remark 2. ./Vote that the condition on t Q from the second part of Theorem^ is really a condition on how 
close the parameters (A*, 7*) are to the "diagonal" given by X*p — 7*0. On the diagonal one trivially sees 
the condition is satisfied. Some algebra shows that the condition is satis fired provided (A*, 7*) lie within the 



v + a \ 'o+i- 1 ■7*0 

l-iL-lfr+q-l-max{p-l,q-l}]) < < 

2pq J X*p 



1 - ^T" ~\P + 1 - 1 - max{p -l,o- 1}] 
2pq 

Theorem 5. Let f(u) = (u + l) p , g(v) := (v + l) q with p, q > 1. 
Suppose 



iV<min<!4 1 " 1 P A ' Q 



p — 1 y p — 1 q — 1 y g- 1 
and (A*, 7*) G T. Then the associated extremal solution o/(JJ)a*, 7 * is bounded. 

Remark 3. Note that p 1— > 4 + + 2y^^y is decreasing and goes to 6 as p — » 00. Hence for N < 6 we 

see all extremal solutions are bounded for any p,q. As in the case of (G)a, 7 one can obtain better results 
provided they restrict the range of the parameters (A, 7) to a certain cone with axis given by J = ^, we omit 
the details. 

We begin with some pointwise comparison results. 
Lemma 3. Let f(u) = (u + l) p and g(v) = (v + l) q where p,q > 1. 

1. Suppose that (u,v) is a smooth solution of (G)\, 1 where Xp > 70. Then v < u < 

2. Suppose (it, v) is the smooth minimal solution of (H)\^ where qX > jp. Then pju > gAi>. 
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Proof: (1) Subtracting two equations of (G)\ n we get 

-A(u-v) = (l+uf-^l + vy-^Xpil + v) - 7 g(l + u)) 

multiply both sides of the above with (u — u)_ to get J \V(u — v)~\ 2 < and therefore v < u. Now, multiply 
the second equation of (G)a. 7 with -2 and again subtract two equations to get 

-A(u - —v) = Ap(l + u) p_1 (l + - «) < 

79 

From maximum principle we get u < ~^v. 

(2) Set K(x) := (u + l) p_1 (t; + l) 9-1 . First note that 

L(u — w) := — A(u — u) — "fpK(x)(u — v) = K(x)((Xq — jp)u + Xq — 7 p), 

and note that the right hand side is nonnegative. If we can show that L satisfies the maximum principle 
then we'd have u — v > 0. 

We now assume that (u,v) is the smooth minimal solution of (H)\^ and additional we assume that 
(A, 7) <E U\T. By Theorem A there is some 77 > and ip > such that 

— Aip — qpyK(x)ip > nip. (15) 

Since (A, 7) ^ T one can infact show that 77 > 0. Hence the linear operator on the left satisfies the maximum 
principle. Since q > 1 we see that L must also satisfy the maximum principle and hence u > v. In the case 
where (A, 7) eWnT we pass to the limit along the fixed parameter ray through (0,0) and (A, 7) and use 
the above result. Hence we have shown that u > v. Now set t := — and then note that 

— 7p 

-A(u - tv) = K(x)Xq((u + 1) - (v + 1)) > 

and hence u > tv. 

□ 

Before we prove Theorem [4] we need a general energy estimate. 

Lemma 4. Let f(u) — (u + l) p and g(v) = (v + l) q . Suppose that (u,v) is a semi-stable solution of (G)a, 7 
with s,t £ R\{±}. Then 

+2pq /(l + uY+P-^l + tj) 5 ^" 1 + p{p - 1) /(l + u) p - 2 (l + vf + q(q - 1) /(l + u) p (l + «) 9 " 2 



<2p(p-l) / (l + u) t+p - 2 (l + u) 9 + 2 M / (l + uj'+f-^l + u)'- 1 
+2 g ( g - 1) J (1 + uf{l + v) s+ «- 2 + 2pq J (I + u) p -\l + vY+i- 1 
Proof. This is an application of LemmaQ] Take <p := (1 + uf — 1 and ip := (1 + v) s — 1 in ([2]), then we have 
pip - 1) /" (1 + U ) p - 2 (1 + w) 9 ((1 + «)* - l) 2 + g(g - 1) /" (1 + u)«- 2 (l + u)* ((1 + w) s - l) 2 
+2pq y (1 + + v)"- 1 ((1 + u)' - 1) ((1 + v) s - 1) 

< - J \Vu\ 2 (l + u) 2t - 2 + J \Vv\ 2 il + v) 2s - 2 (16) 
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Multiply the first and the second equation of (G)\ n with (1 + u) 2t 1 — 1 and (1 + v) 2s — 1, respectively, 
to get 

(2t - 1) / |Vu| 2 (l + u) 2t - 2 = Xp f(l + u) 2t+p - 2 (l + v) q - Xp f(l + uf-^l + v) q 



and 

(2s - 1) / \Vv\ 2 {l + v) 2s - 2 = jq I (1 + v) 2s+q - 2 (l + u) p - jq I (1 + v)^ 1 (1 + uf . 



Using these identities and (|16l) finishes the proof. 



□ 



Proof of Theorem [4j (1) Let (A*, 7*) € T and let (it, v) denote a smooth minimal solution on the ray 
where a := 

Let 1 < t < t + {p — 1) and 1 < s < t+(q — 1) in Lemma 2] to arrive at an inequality of the form 

(U + l) 2t+p - 2 (v + l)' + /(« + 1)> + l) 2s +9- 2 < C t , s . 



First note that 

J |Vu| 2 < pX* J(u + l) p (v + If < C M , 

provided p < 2i+ (p — 1) +p — 2 but this holds since p > 1 and by passing to the limit we see that m* € H&(Q). 
We similarly show that v* g TJo(fi). 

Without loss of generality assume that p > q and hence t+(q — 1) < t+{p~ 1) and so we have 

{u + l) 2t+p+q - 2 < C t , 



for all 1 < t < t+(p — 1). We now re- write the equation as 



A* 

where 



C{ X )U* +p( V * +!)", 



< C(X) = p^ >-{v* + 1)9 < C(u* + l)P+9"2. 

We now apply regularity theory to see that u* is bounded provided c(x), (v* + l) q £ L T for some T > ^. 
But this holds provided 

N 

(p + q-2)- <2t+(p-l)+p + q-2 1 

which is the desired result. To see v* is bounded we now use the pointwise comparison between u and v and 
pass to the limit along the ray T a . 

(2) In (1) we only used the first two integrals from Lemma 2] to obtain estimates. In this part we also 
use the third integral. Let (A*, 7*) € U and let (u,v) denote the a minimal solution on the ray T a , where 
a := p-. The exact proof depends on the sign of X*p — 7*9 and we suppose that X*p > j*q. Let t a < t < T 

and so I p , q ,\* .7*(i), Iq,p,\* ,7*(*) > and p — 1 — ^tzt '9 — 1 — < 0. We now set s = t and examine 
the estimate from Lemma |4j Note the coefficients in front of the first two integrals are negative and the 
coefficient in front of the third integral is positive. The other integrals on the left are lower order terms 
which we drop. Now note that u > v and so we can replace, since the coefficients are negative, the it's in 
the first two integrals from the estimate in Lemma 2] with w's. In the third integral we use the fact that 
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j^(u + I) < v + I. Writing this all out and then again using the fact that we can compare u and v, one can 
see that the following is a lower bound for the left-hand side of the integral estimate given by Lemma [4j 

(K p - 1 - ^) + • (« - 1 - i&i) + 2 « /c + "' 2i+ ' + *" 

+ (, + „ (, + , - 1 - + ((^)' + '-' - l)) /(! + . 

Combining everything gives an estimate of the form 

W-,r(*) /(i + «) 2t+p+? - 2 <c p ,,,a., 7 . ({i + u) t+ 'p+*- 2 . 



Since Ip,q,\* ,y* (t) > we have an estimate. We now proceed exactly as in the first part. We rewrite the 
equation in the alternate form and we then require that 

(p + q-2) — <2t + p + q-2, 

for some to < t where Ip,q,x* ,*/* (t) > 0. 

□ 

Lemma 5. Lei (A*, 7*) G T. 

and let (u, v) denote a minimal solution of (H)\ r( on the ray r CT where a — ¥7. T/ien /or 1 < t < t+(p) 
and 1 < t < t + (q) we have 



1. 



2. 



3. 



where C is uniform on the ray T a . 

Proof. Set (j) := (u + 1)* — 1 and ip := (v + l) r — 1 and put these into the stability inequality given by ([3]) to 
arrive at an inequality of the form 

q( P - ^y) f(u + lf^-Hv + ly- 1 + P ( q j (u + iy- x (v + i)^ 2 - 1 



[(u + lft+P-^v + l)"- 1 <C, 
(u + l)"- 1 (v + l) 2r+q - 1 < C, 
(u + lf+'-^ + l)^- 1 <C, 



+2v/p(p-l)g(g-l) / (« + If +t -\v + l)^- 1 

< C{ P ,q) [(u + l) p+t - 1 {v + l) q - 1 

+c( P ,q) [(u + iyp-^v + iy+T- 1 



and note that for the given choices of t, r the coefficients on the left are positive. One now easily sees that 
the terms on the right are lower order terms and hence we obtain the desired estimates after some standard 
calculations. 

□ 
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Proof of Theorem Without loss of generality we suppose that X*q > ^*p. Let (u,v) denote a 
minimal solution on the ray T a where a := j?. Note that we have u > > v. We first show that 

u* G Hq. First note that 

r |Vu| 2 = Xq J(u + l) p u(v + l)"- 1 , 
along the ray and the right hand side is uniformly bounded provided 

p+Kp-l + 2t+(p), 

which is the case, for any p > 1 and dimension TV, after considering the estimates from Lemma [5] We now 
rewrite the equation for u* as 



-Au* = Xq 



( K+ ^ )P - ) (v* + i)«-V + x q (v* + 



and to show u* is bounded it is sufficient to show that (u* + l) p 1 (v* + l) q 1 g U for some r > ^. Using 
Lemma [5] one sees this is the case provided 

N N 

y(p-l) <p-l + t+(p), -(q-l)<q-l + t + (q). 
So we have shown that u* is bounded and we now use the fact that u* > v* to see the same for v* . 

□ 
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